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ABSTRACT 

A-priori  boundedness  results  for  solutions  of  strongly  coupled  semilinear 
parabolic  systems  of  second  order  under  homogeneous  linear  boundary  conditions 
are  established.  In  contrast  to  [1],  [2],  [4]  it  is  not  supposed  that  the 
diffusion  operator  is  self-adjoint  or  that  the  nonlinearity  is  of  gradient- 
type. 
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SIGNIFICANCE  AND  EXPLANATION 


'Rie  prototype  parabolic  partial  differential  equation  is  the  heat 
conduction  equation.  This  paper  deals  with  systems  of  parabolic  equations. 
Such  systems  occur  in  many  contexts  in  addition  to  heat  conduction,  e.g. 
biology,  in  nuclear  reactor  techniques,  in  economics,  etc. 

Let  the  n-vector  u  denote  the  (unknown)  solution  of  a  system  of  n 
parabolic  partial  differential  equations.  An  Important  question  in  the  study 
of  these  systems  is  the  boundedness  of  u.  Many  techniques  and  criteria  have 

been  developed  to  solve  this  problem  if  the  system  is  weakly  coupled,  i.e.  if 

L  k 

the  k'^^  equation  contains  second  order  space  derivatives  of  only  u'^,  the  k*' 
component  of  u.  If  this  is  not  the  case,  the  system  is  said  to  be  strongly 
coupled . 

In  the  present  paper  for  a  broad  class  of  strongly  coupled  parabolic 
systems  pointwlse  boundedness  of  the  solution  u  is  established.  /" . . . . 


POIMTHISB  A-PRIORI  BOOMDS  POR  STRONGLY  COUPLED 
SEMILINBAR  PARABOLIC  SYSTEMS 


R«lnh«rS  Rftilllnqer 

1.  We  consider  strongly  coupled  semlllneer  parabolic  ayetems  of  the  form 


(1) 


Su" 

3t 


r  3  /  111  3u*\  .  p  hi  3u*  .  . 


in 


for  0<t<T  I  x«fl  ,  k-1,2, 

under  homogeneous  linear  boundary  conditions.  Hera  u  •  (u\...iu")i  0  Is  a  bounded 
domain  In  rf*  and  0  <  T  <  It  la  the  purpose  of  this  paper  to  establish  polntwlse 
a-prlorl  bounds  for  solutions  u  of  (1).  The  same  problem  has  also  been  treated  by 
Cosner  [41 i  AllVakos  (1,  <31  and  Amann  t2«  <71*  In  contrast  to  these  authors/  we  do  not 

VI 

suppose  that  f  Is  a  gradient  or  that  the  matrix  (a^j)  satisfies  any  symmetry 
conditions.  Our  main  result  Is  stated  In  section  2.  Section  3  contains  some  examples. 


2.  Let  net/*  be  a  bounded  domain  whose  boundary  3(1  Is  a  (m-1  )-dlmenslonal  C^-«anlfold 
such  that  n  lies  locally  on  one  side  of  3Q.  Let  0  <  T  <  and  set  J  ■  (0,T)  with 

Jq  •  (OiT).  We  wish  to  study  the  system  of  equations  j 

(2)  “  L^'ct/xlu  +  f’'(t/X/U)  for  t  c  J|  x  e  n,  k  -  1,2,. ..,n  , 

where  u  -  (u\...,u*'),  u^  •  3u’'/3t,  fjJ>'n«rf’-»RP  Is  a  given  nwasurable  function 

and 

1 

L^'ct/xlu  =  T  D  (a^{t,x)D  u*)  +  y  a’^(t,x)D  u*  I 

l,j,t  ’  ^  1,1 

with  -  3/3x^.  Summation  In  l,j  Is  from  1  to  m  and  In  k,l  from  1  to  n.  The 
coefficients  of  l’^  are  assumed  to  satisfy  the  following  smoothness  conditions:  For  all 

1,  j,  k,  t  the  partial  derivatives  Dj  a^*  exist  In  Jq  *  0  and  a^^,  D^  al[^,  aj^  are  j 

bounded  continuous  functions  In  Jq  *  fl.  Moreover,  these  functions  are  H&lder  continuous 
of  exponent  u  In  t  for  some  0  <  u  <  1,  uniformly  with  respect  to  x  t  R.  Also,  the 

_  I 


functions  a][j  are  uniformly  continuous  In  x  with  a  modulus  of  continuity  Independent 

of  t  «  Jq#  and  there  Is  a  function  c  i  Jq  ♦  (0,“)  such  that 

<3)  I  a|[*(t,x)q^q^  >  c(t)  J  (q^)*  In  Jq  *  Q  for  all  q^  e  <?“  . 

Vfe  will  show  at  the  end  of  the  paper  how  condition  (3)  can  be  relaxed. 

We  assume  that  u  satisfies  Dlrlchlet  boundary  conditions, 

(4)  u-0onJx3R. 

y 

However,  If  the  coefficients  of  L  for  all  tc  are  Independent  of  t,  then  we  also  admit 
boundary  conditions  of  the  form  (k  •  1,2,...,n) 

(5)  I  a^J  D.uS^  +  (1-«’')u’‘  +  6*'  J  b'^(x)u*  -  0  on  J  *  9ft  , 

i.j.t  t 

Here  S*'  €  C(9ft,{0,1}),  y  »  (Y •  .y"')  denotes  the  outer  normal  to  9ft  and  the  b''*  are 
continuously  differentiable  functions  on  9ft  satisfying 

(6)  I  fi*'  n''  b*'*  n*  >  0  on  9ft  for  all  n  e  rf*  • 

Of  course,  (4)  Is  a  special  case  of  (5). 

Let  X  »  L^(ft,I^),  2  <  p  <  “,  with  I*!  denoting  the  usual  norm.  For  t  € 

•P  “ 

define  the  operator  A(t)  i  D(A)  C  X  ♦  X  by 

A(t)u  3  -(L''(t,‘ ),...,L"(t,‘ ))U  +  dpU  ,  U  €  D(A) 

with  dp  >  0  a  real  constant  and 

D(A)  "  {u  e  W^'’’(ft,F?*)  ,  u  satisfies  (4)  or  resp.  (5)}  . 

Note  that  0(A)  is  Independent  of  t.  If  we  wish  to  stress  the  fact  that  A(t)  and  X 

depend  on  p,  we  will  write  Ap(t)  and  Xp  In  the  sequel. 

D(A)  is  dense  In  X,  the  A(t)  are  closed  operators  and,  provided  dp  Is  chosen 

sufficiently  large,  we  have 

l(A(t)+X)"^l  <  M  (1  +  |X|)"^  for  all  t  €  Re  X  >  0 

P  ^ 

with  a  constant  Mp  Independent  of  t,  X  (1*1  denotes  the  norm  In  L(X),  the  space  of 
bounded  linear  transformation  on  X).  For  proof  see  (5).  Furthermore,  our  assumptions 
Imply  that  there  are  constants  K,  P  >  0  such  that 

l(A(t)  -  A(3))a"’(t)I  <  Klt-s!*'  ,  t,  8,  T  e  Jp 


-2- 


lA(t)A"’(8)l  <  P  ,  t,  a  e  . 

H«nc«f  for  each  t  e  Jg,  A(t)  generates  an  analytic  semigroup  exp(-8A(t))f  s  >  0,  and 


the  fractional  power 

A®{t) 

of  A(t) 

for  any  a 

>  0  can  be  defined 

as  the  Inverse  of 

A  (t) 

1  r»  - 

“  rial  ^0  • 

.sA(t)sa.1^s 

(see  [9]).  Set  X®  - 

D(a®), 

A  -  A{0), 

with  norm 

1  xl  ■  1  A®  xl  for 
a 

X  e  if* .  We  then  hgve 

the  continuous  Imbeddlngs  (6,  (1«61 

(7a)  X®  ♦  X  for  m/p  <  2a*n/q,  q  >  p  > 

P  <I 

(7b)  X®  ♦  c''(n,rf‘)  for  0  <  V  <  2a-m/p  . 

P 

Finally,  denote  by  H(t,8)/  s  <  t  c  Jq  the  linear  evolution  system  generated  by  the 
operators  A(t),  [9] •  Then  there  Is  a  i  >  0  such  that  for  any  0<8<a<1  we  have 

(8)  IA®{t)W(t,a)A"®(a)l  <  C(t-s)®‘®a"*^'^"“'  ,  s  <  t  £ 

with  a  constant  C  -  C(a,S)  Independent  of  s,  t.  If  the  operators  A(t)  are  Independent 
of  t,  then  W(t,a)  «  exp(-(t-a)A)  and  the  estimate  (8)  Is  wBll-)tnown  (9,  (|2].  In  the 
time-dependent  case  a  detailed  proof  of  (8)  Is  given  In  (8). 

Equations  (2)  vith  (4)  or  (5)  can  be  summarized  In  the  abstract  equation 
(9a)  u,.  +  A(t)u  “  F(t,u)  ,  t  «  J 

where 

F(t,u)  ■■  f(t,*,u)  +  dpU  . 

Together  with 

(9b)  u(0)  -  Uq  e  X 

(9)  describes  an  Initial  value  problem  in  X,  let  us  assume  that  f  satisfies  the 
following  condition: 

(10)  There  are  constants  0  <  5,  p  <  1  such  that  for  any  R  >  0  we  have 

|f(t,x,v)  -  f(8,x,w)|  «  C()t-s|*  +  Iv-wl®) 
for  all  t,  8  €  Jg,  X  £  fl,  V,  w  e  rP  with  a  constant  C  -  C(R). 

It  then  follows  from  (9,  Sect.  2.5]  that  the  initial  value  problem  (9)  has  a  local  solution 
u  e  C((0,T],Xp)  n  ( ( 0,T] , Xp) ,  defined  on  some  Interval  0  <  t  <  T,  provided  Ug  £  1^ 


for  some  6  >  m/2p.  Noraover,  wa  hava  u  t  C((0,t],x^)  for  any  0  <  a  <  B  and 
u(t)  e  0(Ap)  for  0  <  t  <  T.  Wa  wlah  to  aatabllsh  a  tlma-dapandent  bound  for  u  In 
By  [9,  p.  57]  this  is  sufflclant  to  prova  a  global  axlstanca  thaoram. 

'Rius,  lat  u  e  C(Jg,Xp)  n  C^(J,Xp)  danota  from  now  on  a  fixad  solution  of  (9)  In  Xp 
obtalnad  by  tha  raathod  of  proof  usad  In  (9,  Sact«  2.5].  Wa  assuaw  that  p  ~  2  In  casa 
m  <  4  and  that  p  >  m/2  otharwlaa. 

For  t  «  Jq  lat 

o'‘‘(t)  -  nup{(I  Ca!^(t,x)l^)’''2  i  x  e  0) 

1 

and  sat 

a<t)  -  8up{  \  v’'a’^(t)w*  s  V,  w  e  B?',  |v|  -  Iw|  -  1} 

k,l 

with  I* I  danotlng  the  Euclidean  norm  In  bP.  Wa  define 

v(t)  -  r /jj  I  (uNt.xjj^dx  ,  t  e  J 

and  Introduce  tha  following  hypotheses: 

(v^)  There  la  a  continuous  function  4  t  x  ^  such  that 

L  I  u''f’'(t,x,u)dx  <  •(t,V)  In  J  . 

k 

(Vj)  The  maximal  solution  y„„  of  the  ordinary  differential  aquation 

^t  -♦(t,y)  >  y(0)  -V(0) 

with  c  given  by  (3)  exists  on  Jq  and  Is  bounded. 

(Vj)  Thera  are  constants  >  0,  r  >  1  such  that 

If(t,x,z)|  <  Cg(1  +  1*1')  for  t<J,  xetl,  xeB?*  . 

Theorem.  Let  u  be  a  solution  of  (9)  In  Xp  as  described  above  and  assume  uq  e  D(Ap). 
Let  hypotheses  (V^),  (V2)  and  (Vj)  with  r  <  1  +  4/m  be  satisfied.  Then,  for  any  8  <  1, 
u(t)  Is  bounded  In  X^  uniformly  on  Jq.  Xn  particular,  u(t)  la  uniformly  bounded  in 
C(n,B"). 


Remarfcs  (1)  Assume  that  f  satisfies  (10)  and  that  Ug  e  for  some  8  >  m/2p.  (This 
holdSf  if  Uji  e  c'*(5,r")  for  some  u  >  0  with  ug  “  0  on  80  for  all  k  with  6*'  «  0 
in  (5).)  Then  (9)  will  have  a  local  solution  u  in  Xp<  As  noted  above,  u(t)  e  0(Ap) 
for  t  >  0.  He  can  thus  use  the  assertion  of  the  theorem  to  prove  that  this  solution  is 
global.  However,  in  this  case  the  initial  condition  on  y  in  (V2)  has  to  be  replaced  by 
y(0)  •  V(0)  *  e,  where  e  >  0  is  arbitrarily  small. 

(il)  Uniform  a-priori  bounds  on  u  in  L‘*(0,rf')  with  f  satisfying  (10)  have  also  bean 
obtained  by  Cosner  [4]  and  Allkakos  (1,  (3}  for  the  time-independent  system  (2)  under 
Dlrichlet  boundary  conditions  and  by  Amann  (2,  (^7)  for  the  boundary  value  problem  (2), 

(5).  (In  the  last  paper,  time-dependent  boundary  conditions  and  higher  order  elliptic 
systems  are  also  considered. )  The  bounds  on  r  used  in  these  papers  are  less  restrictive 
than  the  one  stated  above.  However,  all  these  authors  Impose  severe  structure  conditions 
on  the  system  (2)i  It  is  assumed  that  a][j  •  a^^  •  a^j  in  Jg  «  0  for  all  1,  j,  k,  t 
and  that  f  is  of  the  form  f  •  q+h,  where  g  satisfies  a  linear  growth  condition  and 
h  is  a  gradient,  h  •  9rad^  H(x,u).  No  such  conditions  are  needed  in  the  above  theorem, 
(ill)  As  the  following  proof  shows  it  suffices  to  require  c(t)  >0  on  Jg  in  (3)  in 
case  a(t)  ••  0  on  Jg. 

Proof  of  the  theorem.  He  first  establish  a  time-independent  bound  for  u  on  Jg  in  Xj. 
In  fact,  since  u  r  C^(J,X2),  we  get  for  t  e  J 


Py  partial  integration  the  first  summand  is  equal  to 


i.,k.  “  ■«  V  ’  ■ ”>"  *‘>  ‘ 

-  I.  -  /« ,  I  .  vM 

<  0  -  c(t)  /_  I  Id  u*'!^  -  -c(t)  /„|grad  u|^  , 

k.j  ’ 


where  we  have  used  (3)  and  (6).  Alao 


^  u’'aj^*'  Dj^u*  <  /jj  a(t)  |u|  Igrad  u| 

l#k#i 


<  -j  v(t)  +  c(t)  Igrad  u|  ^  • 


Summing  ud  and  using  (V^)  we  thus  get 


dV/dt  <  Ht.V)  +  j  V  in  J 


But  this  implies  V  <  “^O'  “''ete  y,^x  defined  in  (Vj).  Hence  V  is  bounded 

in  Jq,  i>e>  u  is  bounded  in  X2>  uniformly  on  Jq. 

Since  u  is  a  solution  of  (9)  in  X2>  it  follows  in  particular  that  u  satisfies 
in  Xj  the  integral  equation 

(11)  u(t)  -  W(t,0)Uo  +  /J  W{t,s)F(8,u(s))ds  '  t  '  Jq  • 

By  (Vj)  and  (7a)  we  have 


IF(8,u(8))I  ,  <  C,{1  +  Iu(s)l\_) 


.2 


,2t 


<  C,(1  t  lu(8)lj) 
2  a 


for  any  u  >  m(r-1)/4r  with  constants  C^,  C2  Independent  of  8  e  J.  We  now  make  use  of 
the  inequality 


lul  <  flul'’,  lulj"''  for  u  e  X®,  0  <  o  <  B 
a  ,2  D  2 


IF(8,u(s))I  2  <  Cj(1  +  Iu(s)lg*’"'’’)  for  s  t  J 


where  u  ■  1  -  0/6  and  C  •  C(a,8)  is  a  certain  constant  (see  [9,  (1.55)]).  This  gives 
(12) 

with  8  >  a  arbitrary  and  constant  Cj.  We  choose  8  in  such  a  way  that 
(13)  a  <  8  <  1  and  r(1-v)  <  1  . 

This  is  equivalent  to  ra  <  8  <  1,  i.e.  to  m(r-1)/4  <  8  <  1.  Since  r  <  1+4/ra,  such  a 


I 


ip 


•>s 


m 

K 


r 

r/A 


•\V. 


Jh  •  • 


i] 


tk 


- -  i  -v-Lr-v-' 
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holes  Is  possible 


Define  s(t)  -  lu(t)lg>  Using  [9,  (1.68)]  one  can  show  that  z  e  C(Jg,X2)>  Further, 
for  t  e  J  we  have 

(14)  s(t)  <  +  /g  C^(t-o)'®e"*‘'^~“’(1  +  z’^'’'^’(8))d8 

by  (8)  and  (12)  with  constants  C^,  Cg  Independent  of  t»  8  and  arbitrary  6  <  ?  <  1. 

Here,  use  Is  made  of  the  estimate 


S  — fi 

IX  A  (t)l  <  Const* 


for 


8  <  8,  t  c  J. 


which  can  be  proved  as  [9,  (1.59)). 

But  r(1-v)  <  1  by  (13)  and  hence  (14)  implies  that  z(t)  <  S  In  Jg  for  some 
constant  S.  In  case  m  <  4  this  establishes  the  assertion  of  the  theorem. 


In  case  m  >  4  we  choose  8  close  to  1  and  apply  (7a).  It  follows  that  u(t)  Is 


bounded  In  X  ,  uniformly  on  Jg,  for  any  q  >  2  satisfying  1/q  >  -r  ~  2/m.  We  can 


•q*  - -  —  -0'  .  2 

choose  q  “  q^  >  2r.  Consider  (9)  In  Xp^  with  p^  «  q^/r.  Then  F(s,u(8))  for  s  £  Jg 
Is  tjounded  In  5^^  by  (Vj)  and  using  (11)  It  Is  easy  to  see  that  lu(t)lg  Is  bounded 
on  Jg  for  any  0  <  8  <  1.  nius,  again  by  (7a),  u  Is  bounded  In  X^^  for  any  qj  >  p^ 
satisfying  2  +  m/q2  >  m/p^.  Repeating  this  reasoning  we  get  two  sequences  (q^  ),  (l\i  ) 
with 

q,. 


“  %  +  1 


<  2  - 


i»+1 


where  q^  >  2r,  1/q^  >  —  -  2/m  and  m  >  4.  It  follows  that  after  finitely  many  steps  we 


can  choose  p^  -  p  >  m/2.  This  proves  the  theorem. 


3.  Examples,  o )  Consider  the  system 

2 

u^  »  X4u  -  p  ^4v  -  V  "  I 

(  15) 

v^  »  -p  jAu  +  pAv  +  uv  -  Y  jV 

under  )x>undary  conditions  (5).  Let  X,U,P,,P2,  I^»72  constants  with 

2  7 

X,  U,  Y,,  Y,  >  0  and  4Xii  >  (p  ♦  p  )  .  Since  g(u,v)  -  (-v‘,uv)  Is  not  a  gradient  and 


since  we  do  not  assume  ~  P^i  results  In  [1],  [2]  and  [4]  cannot  be  applied  to 

(15).  To  prove  global  boundedness  for  solutions  of  (15)  one  could  also  try  to  use 

comparison  methods  [3],  [10].  However,  apart  from  the  structure  of  the  nonlinearity  In 

(15)  this  method  will  only  be  available  If  we  additionally  assume  that  p^«p2"0  or 
2 

that  ^^^2  t  (^~u)  >  0.  In  this  case,  possible  candidates  for  Invariant  sets  are  sets 

which  at  each  boundary  point  have  their  normal  vectors  equal  to  left  eigenvectors  of  the 
diffusion  matrix  (see  [3]).  But  on  the  boundary  of  any  such  set  (at  least  for  small  y 
there  are  points  at  which  the  nonlinearity  In  (15)  does  not  point  Inward.  Hence  this 
method  also  cannot  be  applied  to  (15). 

On  the  other  hand.  It  Is  easy  to  see  that  (V^)  Is  satisfied  with  ♦(t,V)  »  -2yv,  y  “ 
mln(Y<|,Y2)*  Malting  use  of  the  theorem  we  thus  see  that  for  any  smooth  Initial  data  (15) 
will  have  a  bounded  solution  existing  for  all  time  provided  the  space  dimension  m  Is 
1,  2  or  3.  Note  that  y^a^  *  ^  t  ♦  “  and  hence  u,v+0  as  t**. 

6)  Let  0  >  (0,L),  L  >  0  and  consider 


u  “U  +3v  +V-U 
t  XX  XX 


(16) 


+  U  -  V 


with  “x  “  ''x  “  °  Since  the  diffusion  matrix  ^  “  (q  symmetric,  the 

hypotheses  In  [1],  [2]  and  [4]  are  not  satisfied.  Comparison  methods  cannot  be  applied  to 
(16),  since  A  cannot  be  brought  into  diagonal  form.  Also,  A  Is  not  positive  definite. 


and  hence  condition  (3)  Is  not  satisfied. 


1 


Define  w  -  —  u,  *  »  v.  This  gives 


w  ”w  +z  +4-1-  9w^ 
t  XX  XX  3 


(  17) 


+  3w  -  r 


with  “x  ”  *x  ”  transformed  boundary  value  problem,  all  assumptions  of 


section  2  are  satisfied.  In  particular,  we  can  choose 


Hence  (17),  and  thus  (16),  lias  a  global  solution  for  any  smooth  Initial  Yaluea* 

The  transformation  used  In  the  last  example  Indicates  a  way  to  get  rid  of  condition 
(3).  It  Is  only  necessary  to  assume  that  there  are  positive  constants  such  that, 

with  u''  replaced  by  w''  “  the  transformed  boundary  value  problem  (2),  (5) 

satisfies  the  assumptions  of  the  theorem*  In  essence,  this  leads  to  the  following  problem: 
Let  A  be  a  real  ( n*n)-matrlx.  Give  sufficient  general  conditions  on  A  under  which 
there  exists  a  diagonal  matrix  0  -  (d.),*..,d^)  with  positive  dj^  such  that  DA  is 
positive  definite. 

For  n  “  2  the  answer  can  be  given  readily,  but  for  higher  dimensions  the  problem 
seems  to  he  rather  difficult.  A  recent  discussion  of  this  question  and  some  results  can  be 
found  in  [7). 
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